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One of the most enigmatic TeV binary systems, LS I +61° 303 exhibits a high degree of flux 
modulation from radio to TeV wavelengths over a single orbit of the binary system, once every 
~26.5 days. The binary system also exhibits a ~ 4.5 year superorbital modulation in radio, X- 
ray and GeV emission which is yet to be seen in TeV gamma rays. LS I +61° 303 has been 
observed by VERITAS in the TeV energy range and multiwavelength partners (optical - GeV). 
The contemporaneous multiwavelength data sets enable searches for correlations between X-ray 
and TeV emission, as well as GeV and TeV emission. These correlation studies can further 
elucidate the astrophysical properties of this binary system. We present the preliminary results 
from the analysis of recent VERITAS observations from the 2014-2015 season in conjunction 
with 5wiff-XRT (0.3-10 keV) and Fermi-LAT (0.3-300 GeV) observations. 
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1. Introduction 

The high-mass X-ray binary (HMXB) LS I -i-61° 303 is one of the most studied binary systems 
in the gamma-ray sky. It is located at a distance of ~ 2 kpc [1] and consists of a massive BO Ve star 
[2] and a compact object. The nature of the compact object, (either a black hole [3, 4] or a neutron 
star [5]) is a subject of active debate because of the poorly constrained mass of the compact object 
and uncertainties in the inclination angle of the system. No evidence for pulsed emission in the 
radio, X-ray, or GeV bands has been detected thus far. The orbital period (~ 26.5 days [6]) of the 
binary system modulates energetic outbursts in radio [7], X-ray [8], GeV [9] and TeV wavelengths 
[10, 11]. It has a highly elliptical orbit {e = 0.54 ±0.03 [12]) around the Be star which is losing 
mass due to its fast rotation, forming an equatorial disc [2]. The zero orbital phase (0o) is defined 
at MJD 43366, the date of the first radio detection of the system [6]. Based on this the apastron 
passage occurs at an orbital phase (j) = 0.77 and periastron passage at 0 = 0.27 according to radial 
velocity studies [12]. 

The system is also found to exhibit a superorbital period of 1667 ± 8 days in radio [6], Ha 
emission lines [13], X-rays [14] and GeV emission [15]. The source flares in X-ray wavelengths 
as seen by the RXTE satellite between (p ~ 0.35 to 0 ~ 0.75 on the superorbital ~ 4.5 years 
timescale. Radio wavelength bursts are also observed, delayed from the X-ray by A0 ~ 0.2 [14]. It 
is suggested that plasma bubbles consisting of energetic particles travel from a region close to the 
binary system to a radio emission region, located about 10 times the distance of separation between 
the members of the binary system. The time of flight for the plasma bubbles may cause delay of 
the radio emission [14]. 

It was initially observed that the GeV flux as seen by Fermi-LAI is also modulated by the 
orbital period with peak emission at 0 ~ 0.4, near the periastron passage [16]. With further mon¬ 
itoring, it came to light that the periodic GeV emission at periastron occurs only during certain 
epochs. The lightcurve from Fermi-LPH data shows a broader peak after periastron and a smaller 
peak just before apastron [9]. Later timing analysis of the GeV lightcurve demonstrated the pres¬ 
ence of periodicity but not with equal power at all times. There is a discontinuity in the periodicity 
of GeV emission at periastron and intermitent low flux just before apastron [17]. It is to be noted 
that this source exhibits different behavior in the low and high energy range. Although modulated 
by the superorbital period, the radio outbursts occurs consistently at every apastron passage. 

LS 1-1-61° 303 was first detected in the TeV energy range in 2006-2007 by MAGIC [11] and 
confirmed by VLRITAS [10, 18]. Initially, the source was detected in TeV over multiple orbits, 
with periodic outburts near apastron. Between 2008 and 2010 the source was not detected near 
the apastron passage but it was detected at a lower flux level during a single orbit near periastron 
passage [19]. Lurther monitoring in 2011-2012 season with VLRITAS demonstrated LS I -i-61° 
303 was active and flaring in the previously seen apastron region of the orbit [16]. 

The aim of this paper is to present additional multiwavelength monitoring of the source in 
X-ray, GeV and TeV wavelengths by reporting the recent data taken in the 2014-2015 season. 

2. VERITAS Observations 

The Very Energetic Radiation Imaging Telescope Array System consists of four 12m tele- 
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scopes, each with a Davies-Cotton tessellated mirror structure, using the imaging atmospheric- 
Cherenkov technique to observe gamma rays from 85 GeV to > 30TeV. It is located at the Fred 
Lawrence Whipple Observatory (FLWO) in southern Arizona (31°40' N, 110°57' W, 1.3 km a.s.L). 
For detailed description and characterization of the instrument refer to [20, 21]. 



MJD 


Figure 1: Daily light curves for LS I +61° 303 with VERITAS (>300 GeV), Fenni-LAT (0.3-300 GeV) 
and Swift-XRT (2-10 keV) during October 2014 - January 2015. The orbital phase (j) of the system is also 
indicated. Flux upper limits for VERITAS are calculated at the 99% conhdence level for data points with 
significance < 3(J and are marked as downward arrows. Flux upper limits for Fenni-LAIY are calculated at 
the 90% conhdence level and are marked similarly. 

VERITAS has been monitoring LS I +61° 303 every year since 2007. The data presented here 
were recorded from October 16 2014 to January 10 2015 (MJD 56946 to MJD 57032) amounting 
to a total of ~ 28 hours of quality-selected livetime observations. The data is unevenly sampled 
around the apastron passage between (j) ~ 0.4 to 0 ~ 0.9. During the first two orbital phases 
observed in October and November the source was in an extremely active state producing a flux 
of 15 — 70 X 10^^^ ys cm^^ s^^ above 300 GeV corresponding to 25% of the Crab Nebula flux in 
the same energy range. The lightcurve of this source can be seen in Figure 1 . This was historically 
the highest flux recorded ever for this binary system in VERITAS. An excess of 563 events over 
the background was detected corresponding to a statistical significance of ~ 19a. The source was 
found to revert to its typical flux level (~ 15% of the Crab Nebula flux) during the last observed 
orbit in December 2014. As shown in Figure 2, the differential energy spectrum in the range 
0.15-10 TeV is well fitted to a power-law model {x^/n.d.f. = 10.98/6) represented by (2.08 ± 
0.13ifa() X 10^^^ X ( 1 The energy spectrum is found to be 

consistent with previous observations at a lower flux level [10, 16, 19]. Nightly variabilty is also 
observed at the level of 3a level as described in [22] following the method in [16]. 

3. Fermi-LAT Observations 

The Fermi Gamma-ray Space Telescope [23] has two components, the Gamma-ray Burst Mon- 
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Figure 2: VERITAS Spectral Energy 
Distribution > 300GeV during 2014-2015 
observation season. The spectral in¬ 
dex (—2.58 ± Q-QIstat) is consistent with 
previous published value of —2.59 ± 
0.15,, [16]. 



Figure 3: Daily flux correlation study for VERITAS 
and Fermi-LAT. The Pearson correlation coefficient of 
0.071+035 indicates independent data sets. It is con¬ 
sistent with previously published correlation studies. 


itor (GBM) dedicated to the observation of GRBs and the Large Area Telescope (LAT) which has 
a large field of view for gamma rays in the energy range of 100 MeV - 500 GeV. The analysis was 
performed on all data available between October 3rd 2014 (MJD 56933) and January 18th 2015 
(MJD 57040) using the Fermi Science tools vl0r0p5 available from the Fermi Science Support 
Center (FSSC).^ Only "SOURCE" class events (event class 128) were analyzed from Pass 8 data 
with event reconstruction following the standard quality cuts prescribed by FSSC. Other standard 
cuts applied are for zenith angle < 90° to reduce gamma rays from the Earth’s limb. A region of 
interest (ROI) of 15° was initially chosen and all 3EGE sources within it were taken into account 
in the fit. The spectral parameters of all sources beyond 5° from ES I -i-61° 303 were fixed and 
any source wifh a TS < 2 was removed from fhe model. This model was used fo produce daily 
lighfcurves by fixing all sources excepf ES I 61° -i- 303 and fhe nearby pulsar 3EGE J0248.1-1-6021. 

Specfral analysis was performed using likeSED (a user confribufed fool) from ESCC.^ The 
specfrum of fhe source is filled wifh a power law wifh exponential culoff of fhe form 

(3.1) 

IMeV 

where Nq = (9.5 ±3.1) x 10^^ 75 MeV^^ cm^^ s^^, E = 1.97 ±0.05, and ^cutoff = 4.00±0.69 
GeV. 

The correlalion sludy for GeV and TeV flux shown in Eigure 3, yielded a Pearson correlation 
coefficienf of 0.07 l+Qjg indicating independenf dafa sefs as previously observed [16]. 

4. Swift-XKT Observations 

Swift [24] is a spaced-based felescope primarily involved in observing GRBs. If carries Ihree 

'http://fermi.gsfc.nasa.gov/ssc/data/analysis/software/ 
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Swift-XRT PC spectrum of 2CG135+01 
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Figure 4: Swift X-ray spectrum of LS I h- 61° 303 fits a power law model with photoelectric absorption. 


payloads of which the X-ray Telescope (XRT) is capable of imaging and obtaining spectra for 
pointed observations in the energy range of 0.3-10 keV. It has a flexible Target of Opportunity 
program which supports multiwavelength observations. VERITAS made several simultaneous ob¬ 
servations with Stvi/t-XRT for LS I -i-61° 303 during the 2014-2015 season. The Stvtjti-XRT data 
has been analysed using the webtools available from UK Swift Science Data Centre^ . The data is 
processed to create X-ray light curves and spectra based on HEAsoft v6.16 package. 



Figure 5: Correlation study of contemporaneous Swift-XKT and VERITAS flux data points. A Pearson 
correlation coefficient of O.SO^q jg indicates a correlation between the X-ray and VHE gamma-ray fluxes. 
Archival data, shown as open squares showed no evidence of correlation with a correlation coefficient of 
0.36±0.32 [16]. 

The XRT data was recorded between October 4 2014 and January 17 2015 (MJD 56934 - 
MJD 57039). The source spectrum for the complete data set is shown in Figure 4. The best 
fit simple power-law model with photoelectric absorption gives a photon index of 1 .54 ]']q with 
X^/n.d.f. = 611.59/700 and neutral hydrogen density of 0.87]']QQg x 10^^ cm^^. The photon index 
was found to vary from 1.3lQg to ES^qj with a reduced-;^^ value ranging from 0.71 to 1.15 for 

^http://www.swift.ac.uk/user_objects/ 
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individual observations. VERITAS recorded 12 contemporaneous observations with Swift-XRT 
which are marked on the Figure 5. Archival data [16] from 2011/2012 and previous VERITAS 
and MAGIC measurements [25] are also shown in this figure. A Pearson correlation coefficient 
of 0.80 ]]q between Swift-XRI and VERITAS fluxes was found for the current 2014-2015 data. 
The total correlation coefficient considering the current data and archival data was found to be 

U.O/014 . 

5. Results and Discussion 

The results for the multiwavelength monitoring of ES 1-1-61° 303 have been presented for 2014- 
2015. The system was exceptionally bright in TeV during the first two observed orbital cycles in 
October and November 2014 and then retreated to its quiescent flux levels during December 2014. 
A correlation study performed with GeV and TeV gamma-ray fluxes confirmed previous results 
indicating a lack of correlation. On the other hand, an indication of a correlation is seen between 
the X-ray and TeV data sets. No significant variation of the spectral indices has been found in 
X-ray, GeV or TeV gamma rays. 



Figure 6: VERITAS and Fenni-LKT Spectral Energy Distribution 

A detailed GeV-TeV spectral energy distribution (SED), seen in Fig 6 is obtained from Fermi- 
EAT and VERITAS data for 2014-2015 season. Similar to published observations [16], this con¬ 
temporaneous SED demonstrates a cutoff at ~ 4 GeV. No emission is detected between 10-300 
GeV. Below 10 GeV, the data is fitted with a power-law with exponential cutoff. A simple power- 
law is fitted to the data above 300 GeV. The gap in detected emission may be indicative of GeV 
and TeV gamma rays originating from different particle populations. This hypothesis of different 
parent particle population is also supported by the lack of correlation in GeV and TeV fluxes shown 
in Fig 3. 
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A potential correlation seen between the X-ray and TeV gamma ray data suggest that emission 
in these two energies may be produced by the same population of accelerated particles. Leptonic 
models may be a viable option in such a case as suggested previously [19]. Synchrotron radiation 
from high energy electrons may drive the X-ray emission and these electrons may later undergo 
Inverse Compton scattering, upscattering the photons of the companion star to very high energies. 
A similar stronger correlation was seen during a flaring episode in the September of 2007 between 
VHE data from MAGIC and X-ray data from XMM-Newton and Swift [26]. It was postulated 
that uncertainties make it difficult to draw conclusions about correlations outside outbursts unless 
simultaneous observations were recorded. The historically brightest flare of LS I -i-61° 303 seen in 
the 2014-2015 season by VERITAS bolsters the previously found correlation between X-rays and 
TeV gamma rays during September 2007 flaring episode. 
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